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Although CP violation in the B meson system has been well established by the B factories, there 
has been no direct observation of time reversal violation. The decays of entangled neutral B mesons 
into definite flavor states (S" or B'^), and J/ipKl or ccK^ final states (referred to as B+ or -B-), 
allow comparisons between the probabilities of four pairs of T-conjugated transitions, for example, 
B° B- and B_ ^ B°, as a function of the time difference between the two B decays. Using 
468 million BB pairs produced in T(4S') decays collected by the BABAR detector at SLAC, we 
measure T-violating parameters in the time evolution of neutral B mesons, yielding AS+- = -1.37± 
0.14 (stat.)±0.06 (syst.) and AS^ — 1.17±0.18 (stat.)±0.11 (syst.). These nonzero results represent 
the first direct observation of T violation in the B meson system, through the exchange of initial 
and final states in transitions that can only be connected by a T-symmetry transformation. 

PACS numbers: 13.25. Ft, ll.30.Er, 12.15.Ff, 14.40.Lb 



The observations of CP-symmetry breaking, first in 
neutral K decays [ij and more recently in B mesons 
are consistent with the standard model (SM) mechanism 
of the three-family Cabibbo-Kobayashi-Maskawa (CKM) 
quark-mixing matrix being the dominant source of CP 
violation Local Lorentz invariant quantum field the- 
ories imply CPT invariance Q, in accordance with all 
experimental evidence @. Hence, it is expected that the 
CP- violating weak interaction also violates time reversal 
invariance. 

To date, the only evidence related to T violation has 
been found in the neutral K system, where a difference 
between the probabilities of K° — >• if ° and K° K° 
transitions for a given elapsed time (flavor oscillation 
asymmetry) has been measured This asymmetry is 
both CP- and T-violating, independent of time, and re- 
quires a nonzero decay width difference ATk between the 



neutral K mass eigenstates to be observed The 
dependence with AF^ has aroused controversy (ol-fl^. 
In the neutral B and Bg systems, where ATd and AF^ 
are negligible and significantly smaller, respectively, the 
flavor oscillation asymmetry is much more difficult to de- 
tect ■ Experiments that could provide direct evidence 
supporting Tnon-invariance, without being indirectly in- 
ferred from an observation of CP violation, involve either 
nonvanishing expectation values of 7~-odd observables, or 
the exchange of initial and final states in the time evo- 
lution for transition processes. Among the former, there 
exist upper limits for electric dipole moments of the neu- 
tron and the electron The latter, requiring neutrinos 
or unstable particles, are particularly difhcult to imple- 
ment. 

In this letter, we report the direct observation of T vi- 
olation in the B meson system, through the exchange of 
initial and final states in transitions that can only be con- 
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nected by a 7~-symmetry transformation. The method is 
described in Ref. [3], based on the concepts proposed in 
Ref. [11] and further discussed in Refs. We 
use a data sample of 426 fb~^ of integrated luminosity 
at the T{4:S) resonance, corresponding to 468 x 10^ BB 
pairs, and 45 fb^^ at a center-of-mass (cm.) energy 
40 MeV below the T(45), recorded by the BaBAR detec- 
tor [l^ at the PEP-II asymmetric-energy e+e^ collider 
at SLAC. 

In the decay of the T{4S), the two B mesons are in an 
entangled, antisymmetric state, as required by angular 
momentum conservation for a P-wave particle system. 
This two-body state is usually written in terms of fla- 
vor eigenstates, like i?° and S", but can be expressed 
in terms of any two orthogonal states, such as the B+ 
and B- states introduced in Ref. [l5|. They are defined 
as the neutral B states decaying to J/ipK^ and J/ipK^, 
with — >■ TTTT, respectively. The orthogonality is ful- 
filled when there is only one weak phase contributing to 
the B decay amplitude, and CP violation in neutral kaons 
is neglected. Under these conditions, B+ and i?_ are as- 
sociated to CP-even and CP-odd final states, respectively. 

We select events in which one B candidate is recon- 
structed in a 5+ or i?_ state, and the fiavor of the 
other B is identified, referred to as flavor identification 
(ID). We generically denote reconstructed final states 
that identify the flavor of the B as £- X for S° and i+X 
for B'^. The notation (/i, /2) is used to indicate the flavor 
or CP final states that are reconstructed at corresponding 
times ti and t2, where t2 > ti, i.e., i? — ^ /i is the first 
decay in the event and B ^ f2 is the second decay. For 
later use in Eq. ([T|), we define Ar = t2 — ti > 0. We use 
the word "tag" to refer to the fiavor or CP identification 
of the first B to decay (B — >■ /i). At time ti, the fla- 
vor or CP of the second B to decay is orthogonal to that 
tagged by the first decay. The notation i?2(ii) ^2(^2) 
describes the transition of the B which decays at t2, iden- 
tifying its state at ti (orthogonality condition) and ^2 (de- 
cay). For example, an event reconstructed in the time- 
ordered final states {i^ X , J/ip K'^) identifies the transi- 
tion — > B- for the second B to decay. We compare 
the rate for this transition to its T-reversed B- — s> B'^ 
(exchange of initial and final states) by reconstructing 
the final states {J/ip K'^,£~ X). Any difference in these 
two rates is evidence for T-symmetry violation. There 
are three other independent comparisons that can be 
made between B+ B° (J/iPK^J+X), B° B+ 
{£+X,J/^pKl), and S_ ^ B" (J/^KU+X) transitions 
and their T-conjugates, B° B+ (i^X, J/i^K^), B+ ^■ 
B° {J/^jjK°,£-X), and 5° B„ {£-X, J/^pK°), respec- 
tively. Similarly, four different CP (CPT) comparisons 
can be made, e.g., between the B^ i?_ transition and 
its CP (CPr)-transformed B° B_ (B_ ~> B") 15,]. 

Assuming AF^ — 0, each of the eight transitions has a 



general, time-dependent decay rate ^ ( Ar) given by 

g-r,Ar|^ + sin(ATOrfAT) + cos(AmdAT)},(l) 

where indices a — £^,£~ and /3 = iir°,_ft'" stand for 
£'^X,£''X and ccK'^,J/tpKl final states, respectively, 
and the symbol -f or — indicates whether the decay to 
the flavor final state a occurs before or after the decay to 
the CP final state /?. Here, F^ is the average decay width, 
Amrf is the mass difference between the neutral B mass 
eigenstates, and ^ and ^ are model independent 
coefficients. The sine term, expected to be large in the 
SM, results from the interference between direct decay of 
the neutral B to the J/ipK^ final state and decay after 
B^ — B^ oscillation, while the cosine term arises from 
the interference between decay amplitudes with different 
weak and strong phases, and is expected to be negligible 
compared to the sine term. T violation would manifest 
itself through differences between the S"^^ or p values 
for T-conjugated processes, for example between S'i ^0 
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FIG. 1: (color online). Distributions of (a) mES and (b) 
IS.E for the neutral B decays reconstructed in the ccK^ and 
J/tjiK^ final states, respectively, after flavor ID and vertexing 
requirements. In each plot, the shaded region is the estimated 
background contribution. 

In addition to J/i/jK^, B_ states are reconstructed 
through the ■ip{2S)K'^ and XciK'^ final states (denoted 
generically as ccK^), with J/t/j ,'i'{2S) e+e^,^+/i~, 
^{2S) J/iPtt+tt-, xci J/ipl, and K° 

TT+TT^ , tt^tt'^ (the latter only for J/ipK'^). B+ states are 
identified through J/tpKl, with selection criteria identi- 
cal to the most recent CP-violation study by BaBAR [g^] . 
The J/tjjKl candidates are characterized by the differ- 
ence AE between the reconstructed energy of the B and 
the beam energy in the e'^e~ cm. frame, i?boam' while for 
the ccKg modes we use the beam-energy substituted in- 
variant mass mES = Vi^Lam)^ ~ (Pb)^' where p*g is the 
B momentum in the cm. frame. The composition of the 
final sample is determined as in Ref. jioj , through fits to 
the Toes and AE distributions. Figure [1] shows the toes 
and AE data distributions overlaid with the fit projec- 
tions. The final sample contains 7796 ccK'^ events, with 
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purities in the signal region (5.27 < toes < 5.29 GeV/c^) 
ranging between 87% and 96%, and 5813 J/i^K^ events, 
with a purity of 56% in the [AS] < 10 MeV region. 

The flavor ID of the other neutral B meson in the 
event, not associated with the reconstructed i?+ or is 
made on the basis of the charges of isolated primary lep- 
tons, kaons, pions from D* mesons, and high-momentum 
charged particles. These flavor ID inputs are combined 
using a neural network (NN). The output of the NN 
is then divided into six hierarchical, mutually exclusive 
flavor categories of increasing misidentification (misID) 
probability w [2Qi|. Events for which the NN output in- 
dicates very low discriminating power are excluded from 
further analysis. We determine the signed difference of 
proper time Ai = t/j — ta between the two B decays from 
the measured separation of the decay vertices along the 
collision axis. Events are accepted if the reconstructed 
At uncertainty a At and |Af| are lower than 2.5 ps and 
20 ps, respectively. The performances of the flavor ID 
and At reconstruction algorithms are evaluated by us- 
ing a large sample of flavor-specific neutral B decays to 
D(*)-[7r+,p(770)+,ai(1260)+] and J/ipK*°(^ K+tt-) 
final states (referred to as Bflav sample ') i2Ql . The At 
resolution function is the same as in Ref. [20| except that 
all Gaussian offsets and widths are modeled to be pro- 
portional to (TAt- 

We perform a simultaneous, unbinned maximum likeli- 
hood fit to the At distributions for flavor identifled ccK'^ 
and J/ipKl events, split by flavor category. The signal 
probability density function (PDF) is [15| 

Ha,/3(At) oc (7+^(Attruc)i^(Attruc)®7^((5t;aA^)+ (2) 

g-^(-Attruo)i?(-Att,.uo) <E> n{5t- CTAt), 

where Attrue is the signed difference of proper time be- 
tween the two B decays in the limit of perfect At recon- 
struction, H is the Heaviside step function, TZ{5t;aAt) 
with St — At — Attrue is the resolution function, and 
9 a 13 ^""^ given by Eq. (HI). Note that Attrue is equiva- 
lent to Ar (—At) when a true flavor (CP) tag occurs. 
Because of the convolution with the resolution function, 
the distribution for At > contains predominantly true 
flavor-tagged events, with a small contribution from true 
CP-tagged events at low At, and conversely for At < 0. 
Mistakes in the flavor ID algorithm mix correct and incor- 
rect flavor assignments, and dilute the T- violating asym- 
metries by a factor of approximately (1 — 2w). Back- 
grounds are accounted for by adding terms to Eq. ([2]) 
and are incorporated with identical assumptions about 
their At evolution as in Events are assigned signal 
and background probabilities based on the toes or AE 
distributions, for ccifg or J/ipKl events, respectively. 

A total of 27 parameters are varied in the likelihood 
flt: eight pairs of (5'^^,C^^) coefficients for the signal, 
and 11 parameters describing possible CP and T viola- 
tion in the background. All remaining signal and back- 



ground parameters are fixed to values taken from the 
Bfiav sample, J/-0 -candidate sidebands in J/ipK^, world 
averages for and Ato^ [2T|| , or Monte Carlo (MC) sim- 
ulation [i^]. From the 16 signal coefficients we con- 
struct six pairs of independent asymmetry parameters 
(A5±,AC±), (A5±„AC±,), and (A5±,.^, AC±,^), as 
shown in Table U] The 7~-asymmetry parameters have 
the advantage that 7~-symmetry breaking would directly 
manifest itself through any nonzero value of AS^ or 
AC^, or any difference between AS^ and AS'^.^-, 
or between AC^ and AC^j- (analogously for CP- or 
CPT-symmetry breaking). The measured values for the 
asymmetry parameters are reported in Table HI There 
is another two times three pairs of T-, CP-, and CPT- 
asymmetry parameters, but they are not independent 
and can be derived from Table |T] or Ref. 22l |. 



TABLE I: Definition and measured values of the T-, CP-, and 
CPT-asymmetry parameters. These are defined as the differ- 
ences between the (S^ ^„ , ) coefficients, for 5° B- 

and -B+ — > transitions, and those of the corresponding 
symmetry-transformed transitions. The first uncertainty is 
statistical and the second systematic. 



Parameter 



Result 



AS+- 


-^7- 






-1.37 ±0.14 ±0.06 


AS:^ 








1.17 ±0.18 ±0.11 


AC^ 


= ^7- 






0.10 ±0.14 ±0.08 


AC:^ 








0.04 ±0.14 ±0.08 



AS+, = 
AS^ = 
AC+ = 
AC^ = 






-1.30 ±0.11 ±0.07 
1.33 ± 0.12 ± 0.06 
0.07 ±0.09 ±0.03 
0.08 ±0.10 ±0.04 








0.16 ± 0.21 ± 0.09 


AS^j- 


- 9+ 




-0.03 ±0.13 ±0.06 








0.14 ±0.15 ±0.07 


AC^j- 


- r+ 




0.03 ±0.12 ±0.08 








0.55 ±0.09 ±0.06 
-0.66 ±0.06 ±0.04 



c 



t+,K<l 



0.01 ±0.07 ±0.05 
-0.05 ±0.06 ±0.03 



For transition B^ 
asymmetry as 



Ar{At) 



B-, we build the T-violating 
^7-,Ko(At)-H+^^,o(At) 



n 



(At) 



(3) 



where W^JAt) = 'Ha.j3{±At)H{At). With this con- 
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struction, A']-{At) is defined only for positive At val- 
ues. Neglecting reconstruction effects, we have A-r{At) w 

^sm{AmdAt). We build the 



cos{AmdAt) 



other three T- violating asymmetries similarly. Figure [5] 
shows the four time-dependent T-violating asymmetries, 
overlaid with the projection of the best fit results with 
and without T violation [2^ . The results with T invari- 
ance are obtained applying eight restrictions: AS^ = 



AC± = 0, AS"*, = AS-^,^, and AC^ ^ AC^. 



a) 



0.5- 



-0.5- 



-0.5- 



6 8 
At (ps) 



At (ps) 



, c) 
0.5 - +^ 

o*V+^' 



-0.5- 



d) 



-0.5- 




6 8 

At (ps) 



At (ps) 



FIG. 2: (color online). The four independent T-violating 
asymmetries for transition a) i3_ {£'^X,ccKs), b) 

B+ B° {ccK'^,£+X), c) B" B+ {£+ X , J/tP K^l) , d) 
B- — >■ { J/ip K1, X), for combined flavor categories with 
low misID (leptons and kaons), in the signal region (5.27 < 
niES < 5.29 GeV/c^ for ccKI modes and |A£| < 10 MeV for 
J/ipKl). The points with error bars represent the data, the 
red solid and dashed blue curves represent the projections of 
the best fit results with and without T violation, respectively. 



Using large samples of MC simulated data, we deter- 
mine that the asymmetry parameters are unbiased and 
have Gaussian errors. Splitting the data by flavor cate- 
gory or data-taking period give consistent results. Fit- 
ting a single pair of {S, C) coefficients, reversing the sign 
of S under At O -At, or B+ ^ S_ or B° o B" ex- 
changes, and the sign of C under o S'^ exchange, 
we obtain identical results to those obtained in our most 
recent (^-violation study Performing the analysis 

with B decays to ccK^ and J/ipK*^ final states instead 
of the signal ccK° and J/iIjK^, respectively, we find that 
all the asymmetry parameters are consistent with zero. 

In evaluating systematic uncertainties in the asymme- 
try parameters, we follow the same procedure as in (20| . 
with small changes. We considered the statistical un- 
certainties on the flavor misID probabilities. At resolu- 
tion function, and toes parameters. Differences in the 
misID probabilities and At resolution function between 
Bflav and CP final states, uncertainties due to assump- 
tions in the resolution for signal and background compo- 
nents, compositions of the signal and backgrounds, the 
Toes and AE PDFs, and the branching fractions for the 



backgrounds and their CP properties, have also been ac- 
counted for. We also assign a systematic uncertainty cor- 
responding to any deviation of the fit for MC simulated 
asymmetry parameters from their generated MC values, 
adding in quadrature the deviation and its statistical un- 
certainty. Other sources of uncertainty such as our lim- 
ited knowledge of F^, Ato^, and other fixed parameters, 
the interaction region, the detector alignment, and effects 
due to a nonzero AF^ value in the time dependence and 
the normalization of the PDF, are also considered. Treat- 
ing ccKg and J/ip as orthogonal states and neglecting 
CP violation for flavor categories without leptons, has an 
impact well below the statistical uncertainty. The total 
systematic uncertainties are shown in Table |T] [22|] . 




FIG. 3: (color online). The central values (blue point and red 
square) and 2-dimensional CL contours for 1 — CL — 0.317, 
4.55 X 10"^ 2.70 X 10~^ 6.33 x 10"^ 5.73 x 10"^ and 
1.97 X 10~^ , calculated from the change in the value of 2A In C 
compared with its value at maximum, for the pairs of T- 
asymmetry parameters (AS^, AC^) (blue dashed curves) 
and (AS'.^jAC.^) (red solid curves). Systematic uncertain- 
ties are included. The T-invariance point is shown as a plus 
sign (-I-). 



The significance of the T-violation signal is evaluated 
based on the statistical change in log-likelihood with re- 
spect to the maximum (2Aln£). We reduce 2Aln£ by 
a factor 1 -I- maxjmf} = 1.61 to account for system- 
atic errors in the evaluation of the significance. Here, 
mf — 2(ln£ — In >Ci)/s^, where ln>C is the maximum 
log-likelihood, InCi is the log-likelihood with asymmetry 
parameter i fixed to its total systematic variation and 
maximized over all other parameters, and w 1 is the 
statistical change in 21n£ at 68% confidence level (CL) 
for one degree of freedom (d.o.f). Figure [3] shows CL 
contours calculated from the change 2Aln£ in two di- 
mensions for the 7~- asymmetry parameters {AS^, AC^) 
and (AS*.^, AC^). The value of 2 A In C between the best 
fit solution with and without T violation is 226 units 
with eight d.o.f. Assuming Gaussian errors, this corre- 
sponds to a significance equivalent to 14.0 standard de- 
viations (cr), and thus constitutes direct observation of 
T violation. The significance of CP and CPT violation 
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is determined analogously, obtaining 307 and 5 units, re- 
spectively, equivalent to 17a and 0.3(T, consistent with 
CP violation and CPT invariance. 

In summary, we have measured T-violating param- 
eters in the time evolution of neutral B mesons, by 
comparing the probabilities of B'^ — > — s> 
B^ — B+, and i?_ B^ transitions, to their T con- 
jugate. We determine for the main T-violating param- 
eters AS*^ ^ -1.37 ± 0.14 (stat.) ± 0.06 (syst.) and 
AS*.^ = 1.17 ±0.18 (stat.) ±0.11 (syst.), and observe 
directly for the first time a departure from I' invariance 
in the B meson system, with a significance equivalent to 
14(7. Our results are consistent with current CP-violating 
measurements obtained invoking CPT invariance. They 
constitute the first direct observation of T violation in 
any system through the exchange of initial and final 
states in transitions that can only be connected by a 7~- 
symmetry transformation. 
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TABLE I: Measured values of the {S^ ^,0^ coefficients. The first uncertainty is statistical and the second systematic. 







Transition 


Parameter 


Result 


B_ - 


B° 


{j/i<Kl,rx) 




-0.83 ±0.11 ±0.06 


B" - 
B_ - 


^ B- 
B° 


{£-X,ccK°s) 
{J/^PK1,£+X) 




-0.76 ± 0.06 ± 0.04 
0.70 ± 0.19 ± 0.12 


5° - 


^ B- 


{£+X,ccK°) 


Kg 


0.55 ± 0.09 ± 0.06 


B° - 


^B+ 


{rx,j/^Kl) 




0.51 ±0.17 ±0.11 


B+ - 


-> B° 


(ccKlrX) 




0.67 ±0.10 ±0.08 


B° - 


^B+ 


{l+X,J/xpKl) 




-0.69 ±0.11 ±0.04 


B+ ^ 


-^B° 


(ccKle+X) 




-0.66 ±0.06 ±0.04 


B_ - 




{J/^Kl,t-X) 




0.11 ±0.12 ±0.08 


B"- 


^ B_ 


{£-X,ccKl) 


r+ 


0.08 ± 0.06 ± 0.06 


B- - 




{J/^Kl,t+X) 




0.16 ±0.13 ±0.06 




-> B- 


{i+X,ccKl) 


r+ 

Kg 


0.01 ±0.07 ±0.05 


B"- 


^B+ 


{t-X,.]li,Kl) 


W-,KO 


-0.01 ±0.13 ±0.08 


B+ - 


-^B° 


{ccKlrX) 




0.03 ±0.07 ±0.04 




^B+ 


ie+X,J/xbK'i) 




-0.02 ± 0.11 ± 0.08 


B+ ^ 


^B° 


{ccK'^,e+X) 




-0.05 ± 0.06 ± 0.03 



TABLE II: Statistical correlation coefficients for the vector of (5^^, p) measurements given in the same order as in TableU 



100 

1 100 

9 100 

6 100 
18 -14 100 

18 1 -15 100 
27 -16 22 100 

29 -16 21 100 
-37 22 -15 -20 100 

-41 20 -17 -20 100 / 



Only lower off-diagonal terms are written, in %. 
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23 
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TABLE III: Systematic correlation coefficients for tfie vector of (S^ 
Only lower off-diagonal terms are written, in %. 



, iij) measurements given in tlie same order as in Tabled 
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16 


-4 


57 
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100 
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67 
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100 
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10 
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-3 


100 












30 


-19 


57 


59 


10 


58 
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100 










-28 


-10 


39 


13 


43 


21 


-8 


-1 


100 








42 


-20 


60 


68 


57 


72 


-6 


47 


30 


100 






-31 





23 


17 


20 


8 


11 


6 


58 


18 


100 




41 


-27 


70 


66 


46 


64 





71 


32 


81 


20 


100 


31 


-16 


63 


63 


39 


67 


-23 


59 


39 


63 


24 


73 100 


1 


-1 


15 


7 


2 


2 


-31 


5 


23 


7 


5 


18 49 


28 


-23 


73 


72 


52 


61 


-1 


64 


43 


69 


28 


84 83 


V-14 


-13 


12 


-6 


-34 


11 


2 


34 


23 





31 


17 26 



15 15 100 / 




02468 02468 

At (ps) At (ps) 




02468 02468 

At (ps) At (ps) 



FIG. 1: (color online). The four independent CP-violating asymmetries for transition a) — >■ B- {l'^ X , ccK*^) , b) B+ — > 
{ccKl,t+X), c) W B+ {t+X, J/i^Kl), d) B- B° {J/ij} Kl,e.+ X), for combined flavor categories with low misID (leptons 
and kaons), in the signal region (5.27 < mES < 5.29 GeV/c^ for ccKl modes and |AS| < 10 MeV for J/ipKl). The points 
with error bars represent the data, the red solid and dashed blue curves represent the projections of the best fit results with 
and without CP violation, respectively. 
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FIG. 2: (color online). The four independent CPT- violating asymmetries for transition a) B+ — > {ccKs,£^X), b) B+ — >■ 
{ccK°s,£~X), c) B- ^ B° {J/ipKl,£+X), d) B- B° {J/tpKl,rX), for combined flavor categories with low misID (leptons 
and kaons), in the signal region (5.27 < niES < 5.29 GeV/c^ for ccK% modes and |A£;| < 10 MeV for J/ipK°). The points 
with error bars represent the data, the red solid and dashed blue curves represent the projections of the best fit results with 
and without CPT violation, respectively. 




FIG. 3: (color online). The central values (blue point and red square) and 2-dimensional CL contours for 1 — CL = 0.317, 
4.55 X 10"^, 2.70 X 10"^, 6.33 x lO"'', 5.73 x lO"'', and 1.97 x 10"^, calculated from the change in the value of 2Aln£ 
compared with its value at maximum, for the pairs of CP- (left) and CPT- (right) asymmetry parameters [AS^, AC^) and 
{AS^^, AC^j-) (blue dashed curves) and {AS^, AC^), {AS^^, AC^-j-) (red solid curves). Systematic uncertainties are 
included. The CP- and CPT-invariance points are shown as a plus sign (+). 



